Screening induced temperature dependent transport in 2D graphene 
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We calculate the temperature dependent conductivity of graphene in the presence of randomly 
distributed Coulomb impurity charges arising from the temperature dependent screening of the 
Coulomb disorder without any phonons. The purely electronic temperature dependence of our 
theory arises from two independent mechanisms: the explicit temperature dependence of the finite 
temperature dielectric function e{q, T) and the finite temperature energy averaging of the transport 
scattering time. We find that the calculated temperature dependent conductivity is non-monotonic, 
decreasing with temperature at low temperatures, and increasing at high temperatures. We provide 
a critical comparison with the corresponding physics in semiconductor-based parabolic band 2D 
electron gas systems. 



I. INTRODUCTION 



Ever since the successful fabrication of gated two 
dimensional (2D) graphene monolayers and the mea- 
surement of the density-dependent (i.e., gate voltage 
tuned) conductivity of 2D chiral graphene carriers^, 
transport properties of 2D graphene layers have been 
of great interest to both experimentalistaii^'^'^'^'^'iiSiiSi'iS 
and theoristsiiii^ii^iii^ii^. Much of the early inter- 
est focused on the important issue of the scattering 
mechanisms limiting the low-temperature conductivity 
and the associated graphene "minimal conductivity" at 
the charge neutral (Dirac) point. One of the domi- 
nant low-temperature scattering mechanisms^'^ i^'^i"'^'^'^^ in 
graphene is that due to screened Coulomb scattering 
by unintended charged impurities invariably present in 
the graphene environment, e.g. the substrate (and the 
substrate-graphene interface). By reducing the concen- 
tration of charged impurities in suspended graphene^^^, 
very high quality samples (with mobilities exceeding 
200,000 cm^/Vs which is an order of magnitude improve- 
ment over graphene samples fabricated on a substrate) 
have been made. 

There has been substantial recent experimental^iSiiSiiiS 
and theoreticaliiii^iiliMiiSiie Jij^^ig^ ^^^y. ^oth den- 
sity and temperature dependence of graphene carrier 
transport properties. Much of the observed tempera- 
ture dependent graphene properties have been theoret- 
ically studied in the context of the phonon scattering 
mechanisroi^iiiii^iiSii^S which freezes out at low tempera- 
tures. Our theoretical work presented in the current pa- 
per considers temperature dependent graphene transport 
arising entirely out of electronic mechanisms without any 
phonon effects. Our motivation is partially theoretical, 
but we are also motivated by the intriguing recent ex- 
perimental observationiiSiiSiliS of an increasing graphene 
conductivity with increasing temperature at the Dirac 
point, which obviously cannot be explained by phonons 
since phonon scattering necessarily leads to an increas- 
ing carrier resistivity with increasing temperature (as is 
seen in graphene at higher carrier densities away from 
the Dirac point). 



Our theoretical motivation for studying temperature 
dependent graphene transport associated with purely 
electronic mechanisms arises from the extensively stud- 
ied 2D metal-insulator-transition phenomena3ii2^ in 
semiconductor-based 2D semiconductor systems. The 
low-temperature resistivity measurements in conven- 
tional semiconductor based 2D systems (e.g. Si inver- 
sion layers, GaAs heterostructures and quantum wells) 
report the observation of an anomalously strong tem- 
perature dependent (showing an effective metallic be- 
havior) 2D conductivitjsSii^. In conventional 2D struc- 
tures long-range charged impurity scattering dominates 
low-temperature ohmic transport and it is known that 
the temperature dependent screening of charged impurity 
scattering gives rise to the unusual strong temperature 
dependent metallic behavior at low carrier densitiesi^ 
Since in gated graphene layers, which are similar to 2D 
electron systems in confined semiconductor structures, 
the charged impurities are also the main scattering mech- 
anism, one can expect a strong temperature dependent 
conductivity. However, a very weak (less than 5%) tem- 
perature dependent conductivity from IK to room tem- 
perature has been reported for low mobility graphene 
samples J- 

Recently, with more careful measurements in high mo- 
bility samples a strong temperature dependence of car- 
rier conductivity is reportedi^'^'^ We first summarize 
the key experimental features of the measured carrier 
transport in 2D high mobility graphene sample. Experi- 
mentally one finds a density (n*) separating an effective 
metallic behavior (for density n > n*) from an effective 
insulating behavior (n < n*),^- where metal/insulator is 
defined by whether dp/dT > or < 0. The effective 
metallic behavior is characterized by a drop in the tem- 
perature dependent resistivity, p(T), as T — > 0. At low 
density, near the charge neutral Dirac point, the con- 
ductivity of graphene shows a pronounced non-metallic 
T-dependence (i.e., the increase of resistivity with de- 
creasing T). In this paper, we propose a possible theo- 
retical explanation for (at least a part of) the observed 
temperature dependent conductivity in graphene at low 
carrier density. We emphasize that the apparent insulat- 
ing behavior near the Dirac point cannot be explained by 
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phonons. Our explanation is quantitative, microscopic, 
and physically motivated. Our theory is based on an es- 
sential assumption, that is, transport is dominated by 
charged impurity scattering centers (with a density of 
Hi per unit area) which are randomly distributed in the 
graphene environment. We use the finite temperature 
Drude-Boltzmann theory to calculate the Ohmic resis- 
tivity of the graphene electrons. We calculate the 
graphene conductivity in the presence of randomly dis- 
tributed Coulomb impurity charges near the surface with 
the electron-impurity interaction being screened by the 
2D graphene carriers in the random phase approximation 
(RPA). The screened Coulomb scattering is the only scat- 
tering mechanism in our theory. We also compare criti- 
cally the 2D graphene situation with the corresponding 
situation in the conventional parabolic 2D electron gas 
systems. 

We neglect all phonon scattering effects in this calcu- 
lation, which have been considered recently finding that 
acoustic phonon scattering gives rise to a linear resistivity 
with temperaturefi^iii Given that 2D graphene is essen- 
tially a weakly interacting system with effective ~ 0.88 
for graphene on Si02, a constant independent of carrier 
density (rg = e^/nhvp here is just the effective fine struc- 
ture constant), we expect our RPA-Boltzmann theory to 
be a quantitatively and qualitatively accurate description 
of graphene transport for all practical purposes. 

The paper is organized as follows. In section II the 
Boltzmann transport theory is presented to calculate 
temperature dependent 2D graphene conductivity. In 
Section III we study the temperature dependent screen- 
ing function. Section IV presents the results of the cal- 
culations. We conclude in Section V with a discussion. 



by D{e) = g\€\/{2TTh'^Vp), where g = gsQv is the total 
degeneracy {gs = '2,gv — 2 being the spin and valley 
degeneracies, respectively). 

When the external field is weak and the displace- 
ment of the distribution function from the thermal equi- 
librium is small, we may write the distribution func- 
tion to the lowest order in the applied electric field (E) 
fsk = /(esk) + gsk, where /(e^k) is the equihbrium Fermi 
distribution function and gsk is proportional to the field. 
Assuming spatial uniformity and the steady state electric 
field, the Boltzmann transport equation is written as 



d/,k\ dk 9/(e,k) ^ df 



J (27r)^ ~ *^sk,sk' 



where u^k = sui^k/jkl is the electron velocity, and 
Wsk.s'k' is the quantum mechanical scattering probabil- 
ity. Within the Born approximation Wsk,s'k' for scatter- 
ing from s'k' to sk can be written by 

Wsk,s'-k' = |(14k,s'k')l^ ^ (Csk - Es'k') , (4) 

where (^sk,s'k') is the matrix element of the scatter- 
ing potential associated with impurity disorder in the 
graphene environment, and rii the number of impurities 
per unit area. Following the usual approximation scheme 
we have assumed an ensemble averaging over random un- 
correlated impurities. Note that since we consider elastic 
impurity scattering, the interband processes (s ^ s') are 
not permitted. When the relaxation time approximation 
is valid, we have 



II. CONDUCTIVITY IN BOLTZMANN 
THEORY 

The low-energy band Hamiltonian for graphene is well- 
approximated by a two dimensional (2D) Dirac equa- 
tion for massless particles, the so-called Dirac-Weyl 
equation 
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(1) 



where is the 2D Fermi velocity, ax and ay are Pauli 
spinors and k is the momentum relative to the Dirac 
points. The corresponding eigenstates are given by the 
plane wave 



ipsk{r) = ^ exp(ik • r)^;k, 
VA 



(2) 



where A is the area of the system, s = ±1 indicate the 
conduction (-1-1) and valence (—1) bands, respectively, 
and fJ^ = -i=(e*^'',s) with 6'k = tan{ky/kx) being the 
polar angle of the momentum k. The corresponding 
energy of graphene for 2D wave vector k is given by 
Esk = stiF|k|, and the density of states (DOS) is given 



gsk = 7 — eE-v^k- 



de 



sk 



(5) 



where T(esk) is the relaxation time or the transport scat- 
tering time^^, and is given by 

(6) 

where ^kk' is the scattering angle between the scattering 
in- and out- wave vectors k and k'. 

The electrical current density is given by 



J = .9 



d^k 
(27r)2 



evsk/; 



sk- 



(7) 



Using Eq. ([5]) we obtain the conductivity in Boltzmann 
transport theory by averaging over energy 



(8) 



and the corresponding temperature dependent resistivity 
is given by p{T) = l/a{T). Note that /(e^) is the Fermi 
distribution function, /(e^) = {1 -|-exp[(efc — /x)]//cBr}~^ 
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where the finite temperature chemical potential, fJ-(T), is 
determined self-consistently to conserve the total number 
of electrons. At T = 0, /(e) is a step function at the 
Fermi energy Ep = ii{T = 0), and we then recover the 

usual conductivity formula: a = ^-^D{Ef)t{Ef). 

The matrix element of the scattering potential of ran- 
domly distributed screened impurity charge centers in 
graphene is given by 



k,sk' 



e{q) 



1 



(9) 



where q = \k — k'|, 9 = 0kk', and Vi{q) = 2ne'^ /{nq) 
is the Fourier transform of the 2D Coulomb potential 
in an effective background lattice dielectric constant n. 
The factor (1 + cos0)/2 arises from the sublattice sym- 
metry (overlap of wave function). In Eq. ([5]), e{q) = 
£(g, T) is the 2D finite temperature static RPA dielectric 
(screening) function appropriate for graphene^, given by 
e{q,T) = l + Wc((2')n(q,T), where n(g,r) is the graphene 
irreducible finite-temperature polarizability function and 
Vc{q) is the Coulomb interaction. Then, the energy de- 
pendent scattering time T{tk) for our model is given in 
the leading-order theory by 



1 



TTn^ 

h 



(2^ 



Vi{q) 



x(l — cos0)(l -t- COS 6 



S (ek - Ck' 



(10) 



The factor (1 — cos 6*) in Eq. (fTU]) weights the amount of 
backward scattering of the electron by the impurity. The 
(1 — cos 9) factor, associated with the vertex correction 
by impurity interaction in the diagrammatic calculation 
of the conductivity, is always present in transport the- 
ories involving elastic scattering. In normal parabolic 
2D systems the factor (1 — cos 9) obviously favors large 
angle scattering events, in particular, the +kF to —kp 
backward scattering. However, in graphene the large an- 
gle scattering is also suppressed due to the wave function 
overlap factor {l + cos9), which arises from the sublattice 
symmetry peculiar to graphene. The energy dependent 
scattering time in graphene thus gets weighted by an an- 
gular contribution factor of (1 — cos0)(l -I- cos 9), which 
suppresses both small-angle scattering and large-angle 
scattering contributions in the scattering rate. There- 
fore, r is insensitive to both small and large angle scat- 
terings. In fact, the dominant contribution to r comes 
from cos^ 9 — 0, i.e. 9 — 7r/2 scattering, which is equiva- 
lent to the kp "right-angle" scattering in contrast to the 
2kF back scattering in ordinary 2D systems. The impor- 
tance of the right-angle scattering in 2D graphene has 
not been emphasized in the literature. 

We note that there are two independent sources of tem- 
perature dependent resistivity in our calculation. One 
comes from the energy averaging defined in Eq. ([8]), and 
the other is the explicit temperature dependence of the 
finite temperature dielectric function e(q,T) which pro- 
duces a direct temperature dependence through screening 



cos U) 

326. Of 



in Eq. (fTO)) . i.e. r(e) in Eq. (10) also depends exphcitly 
on T due to the dependence of e{q, T) on T. Even if 
T does not have any explicit T-dependence the finite- 
temperature energy averaging of Eq. (8) by itself intro- 
duces a temperature dependence as long as T(e) has some 
energy dependence. For example, if the relaxation time 
T(e) is given by a function of energy e as r oc e", then 
we have a (x T^+". We describe the details of energy 
dependent scattering time and temperature dependent 
scattering time in the following sections. 

Before concluding this basic transport theory section of 
this paper we want to point out the key qualitative sim- 
ilarities and differences in the transport theory between 
2D graphene and 2D semiconductor based parabolic 2D 
systems (e.g. Si MOSFETs, GaAs heterostructures and 
quantum wells, SiGe-based 2D structures) which have 
been studied extensively over the last thirty years22iSi. 
First, the formal Boltzmann theory for carrier transport 
is the same in both systems except for the different angu- 
lar factor, (1 — cos 9) in the conventional 2D systems and 
(1 — cos^ 9) in Eq. (10) for 2D graphene. Formally the two 
theories become identical for isotropic s-wave disorder, 
where the scattering potential is zero-range (i.e. a con- 
stant in the wave vector space)??. But, for the long-range 
Coulomb disorder associated with scattering by random 
charged impurities in the environment, which is of inter- 
est in this work, the factor (1 — cos9) and (1 — 
in Eq. (10) would have very different implication 
course, the explicit differences between 2D graphene and 
2D parabohc systems in the density of states D{e) in Eq. 
(8) and the dielectric function e{q,T) in Eq. (10) would 
lead to different temperature dependent conductivities 
in these two systems even if the angular factors were the 
same. 

The importance of fc^^-scattering in graphene ver- 
sus 2fci?-scattering in the ordinary 2D semiconductor 
systems in determining the transport properties, most 
particularly the temperature dependent conductivity, 
cannot be overemphasized. For example, theoretical 
approaches to understanding the temperature depen- 
dent graphene conductivity using the impurity-induced 
Friedel oscillation a^^i^^ (i.e. the 2kF behavior of the 
polarizability function) immediately run into problem, 
as mentioned above, because unlike the regular 2D sys- 
tems, fci^-scattering, not the 2fci?-scattering, dominates 
graphene transport. In fact, while both approaches lead 
to the prediction of weak temperature-dependent con- 
ductivity in graphene at low temperatures in contrast to 
regular 2D system, the 2kF Friedel oscillation approach 
predicts^^ a weak insulating temperature dependence for 
high-density extrinsic 2D graphene whereas the kp ap- 
proach based on the screening theory used in the current 
work leads to a weak metallic graphene conductivity at 
low T/Tp. This is a qualitative and conceptual differ- 
ence, which applies whenever screening is important in 
determining graphene transport properties. 
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III. TEMPERATURE DEPENDENT 
POLARIZABILITY AND SCREENING 

The important temperature dependence of the scatter- 
ing time r arises from the temperature dependent screen- 
ing in Eq. (10). Thus, before we discuss the temperature 
dependent conductivity we first consider temperature de- 
pendent screening (or dielectric function) i.e., 



e{q,T) = l + v,niq,T), 



(11) 



where n(g, T) is the graphene irreducible finite- 
temperature polarizability function, which is given by the 
bare bubble diagram (calculated at T = in ref. [231 fo^' 
2D graphene) 



n(<?,T) 



4E^^^-^^-'(k'k'), (12) 



kss 



where k' = k-l-q, e^k = shvpl'kl, and /^k = [exp{/3(esk- 
/i)} -|- where the finite temperature chemical poten- 
tial /z(T) is determined by the conservation of the total 
electron density as 



(13) 



where (3 — l/fc^T and Fn{x) is given by 

f°° t"dt 
lo 



Fn{x) 



1 -|- exp(i — x) 
The limiting forms of the function Fi (x) are given by 



(14) 



Fi{x) 



12 



a; In 2 



X 

T 



for \x\ < 1 



(15) 



Fi{x) 



X 

T 



TT 



[x) -I- a;ln(l + e"'^') for \x\ > 1. 



(16) 

Thus we have the chemical potential in both low and high 
temperature limits for graphene as 



H{T) « Ef 



6 \Tf 



41n2 T 



for T/Tf < 1 (17) 
for T/Tf > 1. (18) 



After performing the summation over ss' one can 
rewrite the polarizability as 



where 

n+(<?,T)== - 



n{q,T) = n+{q,T) + n-{q,T), (19) 
/k+ - /k'+](l + COS Okk') 



2L2 



E 



£k - £k' 



+ 



[/k+ /k'+](l - cosOkk') 



Ek + Ek' 



(20) 



and 

n-(g,r) 



—y 

2X2 



[/k--/k'-](l+COsgfcfc. 
Ek - Ek' 



[/k-+/k'-](l-COsgfcfcO 
Ek + Ck' 



(21) 



where Ek — /iw_F|k|, andcos6'fcfc' — (A:-|-(7COS0)/|k-|-q| and 
(j) is an angle between k and q. After performing angular 
integration and using the dimensionless quantities 11 = 
n/Do, where Dq = gEF/inh^vj. is the DOS at Fermi 
level, we have 



n+(g,T) = 



T 



■ln(l 



Ef Tf 

— I dk 

kF Jo 1 exp[/3(efe - /i)J 



VI - (2fc/9)' 



(22) 



and 



n-(g,r) 



8 fcF T> ^ 



F Jo 



1 + exp[/3(efe + ^)] ' 



(23) 



At T = 0, /i(T = 0) = and the Eqs. ^ and ([23 
become the zero temperature polarizabilities, i.e.. 



11+ {q) 
and 



1 _ 



q < 2kF 

iF 

(24) 



n-(g) 



nq 
8k^' 



(25) 



From Eqs. (|22p and (P^)) we have the asymptotic form 
of polarizability at high temperatures (T 3> Tf) 



and at low temperatures [T <ti Tf) 

■2 / rrn N 2 



(26) 



n(g,T) 

n('7,T) 



KT) 
Ef 

JL 

Ef 



l-y(|:) forE,<2M 



+ 



27r2 r2 EF^i 



1 



For g = 2fci? we then have 



Eg . -1 
— sm — 

2^ Eq 

for £q > 2^ 



nq 
(27) 



n((j = 2/c;^,T) 



MT) 
Ef 



TT/i 

2^ 



1 - 



V2 



3/2 



(28) 
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FIG. 1: Temperature dependent graphene polarizability (a) 
as a function of wave vector for different temperatures and 
(b) as a function of temperature for different wave vectors. 



where C,(x) is the Riemann's zeta function. 

To obtain the screening constant or the screening wave 
vector gs, we note that the screened potential 



C/(g) = 



<9) 



nq [1 + v^Xi{q)\ 



(29) 



so that, qs{(i) = qVc{q)^{q) = 2Tre^U{q)/ k. In the q ^ 
long wavelength limit, we then have the finite tempera- 
ture Thomas-Fermi wave vector as 



Qs(T) 



8\n{2)r,kF ( ^ 



for T:^Tf 



6 \Tf 



for T < Tf (30) 



The screening wave vector increases linearly with tem- 
perature at high temperatures (T ^ Tp), but becomes 
a constant with a small quadratic correction at low tem- 
peratures (T < Tp). 

In Fig. [1] we show the finite temperature polarizabil- 
ity n((7, T) (a) for different temperatures as a function 
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FIG. 2; (a) The 2D polarizability function Yl{q,T)/Do, 
which is normalized to the density of states at Fermi level 
(Do ~ gm/h^-K), as a function of dimensionless wave vector 
q/kp, where kp is the Fermi wave vector, for several different 
temperatures T/Tp =0, 0.2, 0.4, 0.6, 0.8, 1.0 (top to bottom). 
In (b) n(g, r)/Do is shown as a function of temperature for 
different wave vectors. The strong temperature-induced sup- 
pression of the 2k p Kohn anomaly in screening is evident in 
the figure even for very low T/Tp. A comparison with Fig. 1 
shows the compelling qualitative difference between graphene 
screening and ordinary 2D screening behaviors. 



of wave vector, and (b) for different wave vectors as a 
function of temperature. Note that for q < 2kp the to- 
tal polarizability has a local minimum near T « 0.45Tp, 
however it increases monotonically for q >2kp. The dif- 
ferent temperature dependence between small wave vec- 
tors {q < 2kp) and large wave vectors [q > 2kp) gives 
rise to very different temperature dependent scattering 
rates for 2D graphene, Eq. (fTU]) . compared to that of 
ordinary 2D systems. In graphene the chiral sublattice 
symmetry suppresses backward (i.e. a scattering induced 
wave vector change by 2kp from -\-kp to —kp) scattering, 
so the temperature dependence of screening at q = 2kp 
is not significant for conductivity while the temperature 
dependence of screening at large-angle scattering {2k p) 
always dominates the temperature dependent conductiv- 
ity in ordinary 2D systems. In graphene we have to 
consider the temperature dependence of polarizability at 



q ^ kp rather than at 5 = 2fci? in order to understand 
temperature dependent conductivity due to the screen- 
ing because kp (rather than 2fci?) scattering dominates 
graphene transport properties. Since the graphene po- 
larizabihty at 2fci? increases monotonically with tempera- 
ture, the temperature dependence of resistivity caused by 
2k p scattering (or equivalently the thermal suppression 
of the Friedel oscillation — the behavior of polarizabil- 
ity at g = 2kp is closely related to Friedel oscillations) 
decreases with increasing temperature. Recently^^ the 
Friedel oscillation was considered in graphene and a lin- 
ear temperature dependent correction to the graphene re- 
sistivity was obtained based purely on the 2fcir-scattering 
in analogy with the corresponding parabolic 2D systems. 
However, as we mentioned above, this correction to the 
resistivity arising from the 2/ci?-scattering is negative in 
graphene in contrast to the regular 2D systems (i.e. in- 
sulating behavior, or the resistivity decreases as temper- 
ature increases), which disagrees with the experimental 
observation. In graphene, since the most dominant scat- 
tering happens at q kp we have to investigate the 
temperature dependent polarizability at q ^ kp, which 
decreases with temperature (for T <C Tp) and gives rise 
to the increases of resistivity with temperature as ob- 
served experimentally. 

Now we compare the temperature dependent polariz- 
ability of graphene with that of ordinary 2D systems. 
In Fig. [2] we show the corresponding parabolic 2D po- 
larizability normalized by the density of states at Fermi 
level, Dq — gm/lT?2Tr, where 17 is a degeneracy factor and 
TO is the effective mass of electron. Note that the tem- 
perature dependence of 2D polarizability at q = 2kp is 
much stronger than that of graphene polarizability. Since 
in normal 2D systems the 2k p scattering event is most 
important for the electrical resistivity, the temperature 
dependence of polarizability at q = 2kp completely dom- 
inates at low temperatures (T <C Tp). It is known that 
the strong temperature dependence of the polarizability 
function at q = 2kp (see Fig. [5]) leads to the anoma- 
lously strong temperature dependent resistivity in ordi- 
nary 2D systems. However, the relatively weak temper- 
ature dependence of graphene polarizability for q ^ kp 
compared with the 2D polarizability function at q — 2kp 
should lead to a weak temperature dependent resistivity 
in graphene for T <C Tp. 

Finally, an analytic comparison between graphene and 
2D parabolic systems completes the comparison between 
2D graphene and 2D parabolic semiconductor screening 
properties, as shown in Figs. 1 and 2 respectively. We 
provide below the low (T < Tp) and the high (T > Tf) 
temperature analytic limits for the regular 2D polariz- 
ability function in both the q — Q Thomas-Fermi and 
the q = 2kp Friedel oscillation regimes (to be contrasted 
with the corresponding graphene formula given in Eqs. 
(26) - (28) above). For T < Tf 



U2Di2kp,T)^ D2D 

and for T :s> Tp 



_^Tp 
6kl T 



(32) 



(33) 



where D2D = gm/2'KlT? is the regular 2D DOS. We note 
that in the T 3> Tp limit, the Fermi surface is completely 
thermally suppressed, and therefore the 2fcF-screening or 
2A;F-Friedel oscillation does not carry any special signifi- 
cance, leading to the same T 3> Tp asymptotic screening 
formula in Eq. (33) for all wave vectors. For q = 0, in 
the T Tp limit, we get the usual Debye screening for 
the regular 2D electron gas system, which follows from 
putting (7 = in Eq. (33): 



H2D(g = 0,T»TF)»i?2D 



Tp 
T 



(34) 



For graphene the corresponding high-temperature 
screening formula can be easily derived to be 



H(9,r»TF)«i?0 7^ 
J p 



ln4 



24/c2 



Tp 
T 



(35) 



A comparison of Eqs. (34) and (35) show that the 
high-temperature Debye screening behavior are differ- 
ent in graphene and regular 2D systems just as the low- 
temperature screening behaviors also are. 

In the next section we calculate the temperature- 
dependent graphene conductivity due to the scattering 
by screened Coulomb impurities using the temperature 
dependent screening properties calculated in this section. 



IV. CONDUCTIVITY RESULTS 

A. (y{T) due to screening 

Using the temperature dependent screening wave vec- 
tor, qs{T) of Eq. (30), we can calculate analytically 
the temperature dependent scattering time of charged 
Coulomb impurities arising purely from screening in sys- 
tematic T /Tp <C 1 or ^ 1 asymptotic expansions in the 
low and high temperature limits: 



1 



t{T) 2Trh {hvpf Jo 

_ n, (2TTe^Y 2 
27rfi, V K 



d0{l - cos^ ( 
2 .1 



dx- 



e{q,Ty 



SkJo [x + qs{T)/2k] 
At low temperatures {T -^Tp) we have 



^36) 



7 



where tq = t{T = 0) is the scattering time at T = and 
given by^® 



and 



1 

To 



_nj_ 

2TTh 



2io 



(38) 



dx 



2+n 



(x + 2rs) 
1 + 71 9rs 

At high temperatures {T ^ Tp) we have 

1 1 TT / 1 



t(T) to 16/o V41n(2)r 



(39) 



(40) 



Then from Eq. ([8]) we have the temperature dependent 
conductivity due to screening at low temperatures (T <C 

Tf) 



3 h \TfJ 



(41) 



where ctq = e'^VpD{EF)TQ/2. The calculated conductiv- 
ity decreases quadratically as the temperature increases 
and shows typical metallic temperature dependence. On 
the other hand, at high temperatures [T/Tf ^ 1) we 
have 



^.l^[41n(2)..]^. . 

(To IT J 



T Y 



(42) 



The temperature dependent conductivity due to screen- 
ing effects increases as the temperature increases in the 
high temperature regime, characteristic of an insulating 
system. We note that the temperature dependence is 
weak for T < Tf and is strong for T > Tp- 



B. o-{T) due to energy averaging 

Let us now discuss the temperature dependent conduc- 
tivity due to energy averaging. For graphene the energy 
dependent scattering time can be expressed by 



1 



2k 



' dq 



1 - 



vdq? 



-(e) 2TTh{hvFf Jo kk^\^ \2kJ e[q,TY' 

(43) 

where = hvph and Vi{q) the impurity scattering poten- 
tial. For the unscreened Coulomb potential e(g, T) = 1 
and Vi{q) = 2Tie^ j nq. Thus we have 



T{e) 



1 Ef 



where 



1 Ui / 27re^ 

Tl 4:h\ K 



1 



(44) 
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FIG. 3: Temperature dependent conductivity for unscreened 
Coulomb potential (top line) and screened potential (bottom 
line) in the low temperature limit. 



Then from Eq. (8) we have 



a(T) 



2T2 



(46) 



where ai — e'^VpD{EF)Ti/2. Using Eqs. (15) and (16) 
we have the conductivities for the unscreened Coulomb 
potential scattering in the low temperature limit (T <ti 
Tf) 



0-1 



7r2 

3 \Tf 
ai [l + 0[{T/TFf]], 



and in the high temperature limit [T ^ Tp) 



(47) 



(48) 



Therefore the temperature dependent conductivity, as 
implied by energy averaging only, is almost a constant 
in the low temperature limit (See the top line in Fig. [3]). 
But at high temperatures the conductivity due to energy 
averaging increases as similar to screening effects. 

Now consider the screened Coulomb potential scatter- 
ing. Expanding Eq. (I43p with respect to (e — Ep) we 
have 



1 



r{e) 



To Ep 



1 — fll 



£k — Ep {sk — EpY 



Ep 
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where ai = 2Ji/Jo, a2 = iJ^/Jo, where 

c2+"Vl - x"^ 



J n 



dx- 



{x + qo) 



2+n 



(49) 
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FIG. 4: Calculated resistivity arising from the short-range 
white-noise disorder scattering as a function of scaled temper- 
ature T/Tf. Here ps{T) = l/as{T) and pso = 1/<Js{T = 0). 



where go — iTp/'^kF = Zr^. Then from Eq. (8) we have 



(^0 



oe 



e-Ep ,2 A^-Epf 
1 + ai — ;^ h (a^ — a2)- 



(51) 



With rs = 0.88 we have ai = 0.52 and 02 = 0.22. There- 
fore, the temperature dependent conductivity due to en- 
ergy averaging becomes 



cr(T) = (To 



f T 



(52) 



For screened Coulomb potential the conductivity shows 
metallic behavior, again the temperature dependent cor- 
rection being quadratically weak in the small parame- 
ter t = T/Tp <C 1. Combining results from IV A and 
B, we conclude that the T/Tp <C 1 graphene conduc- 
tivity will have weak quadratic metallic corrections, and 
for T ^ Tp the conductivity increases proportional to 
{T/Tpf. 

Finally in this subsection, we comment briefly on the 
temperature dependence of the conductivity arising from 
the short-range white-noise disorder scattering, which 
may be important in graphene (as well as regular 2D 
parabolic electron systems in semiconductor heterostruc- 
tures) at high carrier densities, where Coulomb disorder 
effects are typically screened out. The temperature de- 
pendence of conductivity for both graphene and regu- 
lar 2D systems due to unscreened short-range disorder 
arises now entirely from the energy averaging effect by 
definition since the short-range disorder is considered 
unscreened. It is easy to calculate (Ts{T) of graphene 



due to the short-range disorder (with scattering strength 
Vi = vq) by carrying out the appropriate energy averag- 
ing, we have 



1 



(53) 



where CTso = e'^VpD{Ep)Ts/2 with Tg — jf-Epv^/ (hvp)^ . 
The following analytic asymptotic results are obtained: 



(t(T < Tp) w (7gQ 



a{T»Tp)^^ 



1 



-Tf/T 



Tp 



Sln2 V T 



(54) 



(55) 



In the low temperature limit the temperature dependence 
of conductivity is exponentially suppressed, but the high 
temperature limit of the conductivity approaches aso/2 
as T ^ 00, i.e., the resistivity at high temperatures in- 
creases up to a factor of two compared with the low tem- 
perature limit resistivity. In Fig. [?] we show the calcu- 
lated resistivity due to the short range disorder scatter- 
ing. 

For the sake of completeness, we provide below the 
equations describing the asymptotic low^ and high— 
temperature behaviors of 2D conductivity for the usual 
gaped parabolic 2D electron system as found in Si MOS- 
FETs and GaAs heterostructures: 



cr(T < Tp) W (7, 



2D 



T 
TV 



C2 



T 



2/3' 



a{T ■:^Tp)Ka 



2D 




(56) 



(57) 



In Eq. dini), af^" = (t(T = 0), and Ci = 2go/(l + go), 
C2 = 2.65g^/(l + go)^, where go = qTp/^kp {qrp 
and kp are the 2D Thomas- Fermi wave vector and 
Fermi wave vector, respectively). In Eq. ([57)l . af^ = 
{e^ /h){n/ni)TrqQ. We note that in the parabolic 2D sys- 
tem qtp = gsgv/a-B, where as is the Bohr radius. We 
have assumed an ideal 2D electron gas here with zero 
thickness in order to compare with the 2D graphene sheet 
which also has a zero thickness. 

In comparing graphene temperature dependence with 
the regular parabolic 2D system, we note the follow- 
ing similarities and differences: (i) For T <C Tp, both 
graphene and parabolic 2D systems manifest metal- 
lic temperature dependent conductivity; (ii) the low- 
temperature (T ^ Tp) conductivity manifests much 
stronger linear temperature dependence in the parabolic 
2D system compared with the quadratic temperature de- 
pendence in graphene; (iii) at high temperatures (T ^ 
Tp), both systems manifest insulating temperature de- 
pendence, but with different power laws in temperature. 
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FIG. 5: Calculated resistivity as a function of scaled temper- 
ature T/Tf for different = 0.88, 2.2, 0.1, 0.01 (from top 
to bottom). Ta = 0.88 (2.2) corresponds to graphene on the 
Si02 substrate (in vacuum). Inset shows the magnified view 
in the low temperature limit T < 0.5Tf- 



C. Numerical results 

In this section we present our directly numerically cal- 
culated resistivities of Eq. ([8]) incorporating all effects 
discussed in previous sections. Our numerical results 
agree completely with the analytic results given in Sec. 
IV A and B of this paper in the appropriate T/Tp <C 1 
and 3> 1 limits. 

In Fig. [5] we show our calculated resistivity as a func- 
tion of temperature for different values. Here — 0.88 
(2.2) corresponds to graphene on Si02 substrate (in vac- 
uum). The small values of Ts, independent of carrier 
density and representing the fine structure constant of 
graphene indicate a weak-coupling system in terms of 
eletron-electron interaction. Note that the calculated 
p{T)/p{T = 0) scales for all electron densities, and there- 
fore results can be shown as a function of a single dimen- 
sionless temperature variable T/Tp for a specific value, 
i.e. p{T) = p{T/Tp;r,). Thus, Fig. [S]can be applied for 
all graphene samples for a given value. In the low 
temperature limit the calculated p(T) increases weakly 
quadratically with temperature, manifesting metallic be- 
havior. But at high temperature p{T) decreases quadrat- 
ically. Thus, we find that the calculated resistivity shows 
a non-monotonicity, i.e., at low temperatures the resis- 
tivity shows metallic behavior and at high temperatures 
it shows insulating behavior. The non-monotonicity of 
temperature dependent p{T) can be understood from the 
screening behavior: the temperature dependent polar- 
izability of graphene shows non-monotonic behavior for 
q < 2kp. The metallic behavior is the strongest at w 1, 
and it becomes weaker as decreases. For > 1 the 
strength of metallic behavior decreases very slowly (see 
the inset of Fig. 5). 

In Fig. [n]we show the calculated temperature depen- 



FIG. 6: Calculated conductivity for different temperatures 
T = 0, 200, 300, 500K (top to bottom) as a function of density. 
We use rs = 0.88 and an impurity density rii — 5x lO^'^cm"^. 




FIG. 7: p{T)/p{0) of an ordinary 2D system for different 
Ts values as a function of temperature. As increases the 
metallic behavior becomes stronger. 



dent conductivity for different temperatures as a function 
of density. In the high (low) density limit the conduc- 
tivity decreases (increases) as the temperature increases. 
Therefore the conductivity shows a non-monotonic be- 
havior, i.e. (j{T) has a local minimum at a finite temper- 
ature. 

For comparison we show, in Fig. [7l the calculated tem- 
perature dependent resistivity of ordinary 2D systems 
for different interaction parameters (= me'^ / Ky/nn), 
which for parabolic 2D systems, in contrast to graphene, 
depend on carrier density. We have used the temperature 
dependent polarizability of Fig. 2 in this calculation. Un- 
like graphene (Fig. [5]) the scaled temperature dependent 
resistivity of ordinary 2D systems depends strongly on 
the electron density (or rs). Since the most dominant 
scattering occurs at q — 2kp and the temperature de- 
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pendence of screening function at 2kp is strong, the cal- 
culated 2D resistivity shows the strong anomalous linear 
T metallic behavior, which is observed in many different 
semiconductor systems (e.g., Si-MOSFET^^, p-GaAa^i, 
n-GaAs^^, SiGe^^, AlAs^). In contrast to the ordinary 
2D systems the observed resistivity of graphene shows 
very weak temperature dependence in high density and 
low mobility samples^i^i^. It has been reported that the 
measured resistivity change of low mobility high density 
samples is less than 10% between 5K and 30QK if one 
takes out the phonon contribution. The weak tempera- 
ture dependence of graphene resistivity can be explained 
by the weak temperature dependence of the screening 
function. Qualitatively, however, the temperature de- 
pendence of graphene resistivity and that of a regular 
semiconductor-based parabolic 2D system is similar from 
the perspective of a large change in T jTp. The calculated 
/3(r) for a regular parabolic 2D electron gas system in 
the presence of screened Coulomb scattering in the range 
T jTp = — 3 also shows the non-monotonicity appar- 
ent in Fig. [71 albeit at somewhat higher (rj-dependent) 
values of T/Tp. In graphene, see Fig. El the resistivity 
maximum typically occurs around T jTp < 0.5 compared 
with T/Tp ~ 1 — 2 for parabolic 2D systems as shown in 
Fig. El 

Given the great deal of work on the temperature- 
dependent transport properties of semiconductor-based 
parabolic 2D electron systems over the last 15 yearsSSi, 
it may be useful for us to discuss the similarities and 
differences between these two 2D systems (i.e. graphene 
and regular 2D electron gas) with respect to the temper- 
ature dependence of the resistivity arising from screened 
Coulomb scattering. The most important qualitative dif- 
ference is that the leading low-temperature {T/Tp <C 
1) correction to the resistivity p(T) in ordinary 2D 
(graphene) systems is linear (quadratic) in tempera- 
ture — both are metallic corrections, i.e. p(T) ^ 
po [1 + 0{T/Tp)] in ordinary 2D systems and p{T) ^ 
po [l -I- 0(T/Tf)^] in graphene. This important qual- 
itative difference, of course, leads to a huge quantita- 
tive difference between the two systems in the sense that 
graphene manifests much weaker temperature-dependent 
resistivity than 2D semiconductor systems at low T/Tp. 
This difference in the quantitative temperature depen- 
dence of the low-temperature resistivity (i.e. linear 
and strong in 2D systems, and quadratic and weak in 
graphene) is further exacerbated by the fact that the 
Fermi temperature in extrinsic graphene tends to be very 
high (e.g. Tp > lOOOK for n > lO^^cm"^ in graphene) 
compared with that of semiconductor-based parabolic 2D 
systems (e.g. Tp ~ 7K for n ~ lO^cm"^ in Si MOS- 
FETs), leading to much smaller effective values of T/Tp 
in gated graphene in the extrinsic high-density regime. 

The above discussion and the results of Fig. explic- 
itly establish that extrinsic graphene should manifest a 
very weak screening induced metallic temperature depen- 
dence in its low-temperature resistivity. This fact is in 
excellent agreement with the available experiments ex- 




100 200 300 

T(K) 

FIG. 8: p(r)/p(0) for n = 5 X lO" cm"^ and for different 
Dingle temperatures (level broadening) Td = 0, 0.05, 0.1, 
O.ISTf as a function of temperature. 

cept near the Dirac point where the system has very low 
carrier density (and is almost intrinsic in nature). Tak- 
ing into account that rF(graphene)« 1500\/nK, where 
h is the 2D graphene carrier density measured in units 
of lO^^cm"^, it is obvious that the screening contri- 
bution to the temperature dependent conductivity of 
graphene, going as {T/TpY, would be extremely small 
in the T — Q — lOOK regime except in the intrinsic 
regime where n 1. The strong metallic temperature 
dependence of the low-temperature conductivity, which 
has been much discussed in the context of the 2D metal- 
insulator transition phenomenon in parabolic 2D semi- 
conductor systems, is therefore absent in gated graphene. 

We now discuss the implications of our theory for 
graphene transport at (or near) the charge neutral Dirac 
point where the carrier density is very low. Since 
(j{T/Tp) [or p{T/Tp)] is a universal function of T/Tp 
in the screening theory, the only difference between the 
low-density Dirac regime and the high-density extrinsic 
regime arises from the effective value of T/Tp due to 
the facts that Tp oc yjn, and that the T/Tp <^ 1 and 
T/Tp > 1 regimes manifest qualitatively different tem- 
perature dependence (c.f. Fig. [5]). For example, for 
n < lO^^cm"^, Tp < lOOif, and it is entirely possi- 
ble for such a low-density regime to manifest the high- 
temperature 'insulating' temperature dependence, as is 
apparent for T/Tp > 0.5 in Fig. O where p{T) decreases 
with increasing temperature dependence. Close to the 
Dirac point Tp (oc ^/n) is arbitrarily small, and there- 
fore the recent experimental observation^ - of a decreas- 
ing p(T) at low T at the Dirac point may simply be a 
reflection of this high-temperature 'insulating' behavior 
arising in our screening theory. The fact that this ap- 
parent insulating temperature dependence is observedi^ 
only at low densities is further evidence in support of the 
screening scenario. 

One puzzling issue in this context is that our theory 
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would predict a metallic cr(T) for T/Tp ^ 1 indepen- 
dent of carrier density except that the T /Tp < 1 regime 
necessitates going to lower temperatures at lower den- 
sities. This seems not to be experimentally observedlS 
near the low-density Dirac regime where an insulating 
temperature dependence is reported at low gate voltage 
down to the lowest temperatures. One reason for the ab- 
sence of the metallic regime could be the suppression of 
the temperature dependent screening by impurity scat- 
tering induced level broadenings as was originally sug- 
gested by one of us some time ago^. Essentially, for 
T < Td, where To = wT/kB is the so-called Dingle 
temperature associated with the impurity-scattering in- 
duced level broadening F, the temperature dependence 
of screening is suppressed by scattering effects. Such a 
broadening-induced suppression of the temperature de- 
pendence of the screening function Il{q,T) for T ^ Td 
can be theoretically incorporated in our conductivity cal- 
culations, and we discuss below such a scenario. It is clear 
that this broadening-induced suppression of the temper- 
ature dependence of screening will strongly suppress the 
metallic behavior of the conductivity for T < To, and 
may qualitatively explain the experimentally observed in- 
sulating cr(T) near the Dirac poin1^. 



In Fig [5] we show the level broadening effects on 
p(T)/p(0). Without level broadening (i.e.. To = 0) 
the low temperature metallic behavior of resistivity is 
quadratic. However, the quadratic temperature depen- 
dence is cut off at low temperatures due to the round- 
ing of the sharp corner in the 2D screening function 
by impurity scattering effects at very low temperatures 
T < Tu, and the explicit temperature dependence of 
e((7, T) is suppressed. Thus, the temperature dependence 
of graphene resistivity with the level broadening included 
in the screening function becomes effectively linear at low 
temperatures. 



In Fig. [5] we show the temperature dependent resis- 
tivity for different densities up to room temperature. In 
the low density limit the Fermi temperature is low, so we 
can see both metallic (T <C Tp) and insulating {T >Tp) 
behaviors if we neglect collision broadening. However, 
a finite Dingle temperature suppresses the low tempera- 
ture metallic behavior at low densities if To > 0.45rF, 
and the system manifests only the insulating dp/ dT < 
temperature dependence for all temperatures. For high 
density samples (To -^Tp) we see only the metallic be- 
havior in Fig. [51 and the level broadening gives rise to 
the linear behavior of the resistivity instead of quadratic 
behavior at low temperatures. It is possible that this 
scenario is operational in the experiments of ref. (lo| 
where an insulating temperature dependence of resistiv- 
ity is observed in graphene at the low-density intrinsic 
regime where To > Tp may apply. 
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FIG. 9: Temperature dependent resistivity for different den- 
sities, n = 5, 10, 20, 50, 100, 500 x 10^° cm"^ as a function of 
temperature (a) without collisional broadening and (b) with 
collision broadening of Tb = lOO/C'. 



V. CONCLUSION 

We have developed a detailed microscopic transport 
theory for 2D graphene conductivity at finite temper- 
atures, assuming charged impurities as the dominant 
source of scattering and neglecting all other scattering 
sources (e.g. phonons). We find that the temperature 
dependant resistivity induced by the temperature de- 
pendent screening is non-monotonic. It shows metallic 
behavior at low temperatures, but insulating behavior 
at high temperatures. The quadratic temperature de- 
pendent correction to p(T) at low temperatures is sup- 
pressed by level broadening effects which give rise to an 
effective linear temperature dependence of resistivity for 
T <C Tp. Even though our Drude-Boltzmann transport 
theory explains both the high density metallic behav- 
ior and the low density insulating behavior as observed 
experimentally, we emphasize that the theory is strictly 
quantitatively valid only in the relatively high density 
regime where our linear screening theory based on the 
homogeneous carrier density model is valid. In the low- 
density regime, near the Dirac point, it is well estab- 
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lished that the graphene layer becomes spatially inhomo- 
geneous with random charged impurity induced electron- 
hole puddles dominating the carrier density profilei^^. 
In this low-density inhomogeneous regime our linear 
screening theory, based on an average density approxima- 
tion, is at best of qualitative validity. It is, therefore, not 
surprising that, although we obtain a reasonable quan- 
titative agreement with experiment in the high-density 
regime away from the Dirac point where p{T) manifests 
very weak metallic temperature dependence, our results 
are only in qualitative agreement with the experimental 
data in the low-density Dirac point regime where p{T) 
shows an insulating temperature dependence^. Loosely 
speaking, our screening theory is valid for n > so 
that the impurity-induced puddle formation is weak, but 
we believe that the qualitative behavior predicted by 
our theory has rather broad validity. Indeed, our the- 
ory provides a plausible qualitative explanation for the 
observed^ weakly metallic and strongly insulating be- 
haviors of the temperature dependent resistivity at high 
and low densities respectively. 

We conclude by pointing out that we have only con- 
sidered in this work the contribution to the temperature 
dependent graphene conductivity arising from screened 
Coulomb disorder, including both the explicit temper- 
ature dependence of the screening function and the im- 
plicit temperature dependence due to the thermal energy 
averaging in the Boltzmann theory. There are other scat- 
tering mechanisms contributing to the temperature de- 
pendence of carrier transport properties, most notably, 
phonon scattering which we have studied elsewherei^. 
When the temperature dependence of the conductivity 
is weak, Matthiessen's rule should apply giving p{T) = 
Pi{T) + pph{T), where Pi(T), Pph(T) are respectively the 
graphene resistivity due to charged impurity scattering 
and phonon scattering. Since pph{T) is very strongly 
suppressed at low temperatures due to the well-known 
Bloch-Griineisen behavior, it is reasonable to expect that 
p{T) is dominated by Pi{T), considered in this work, 
for T < 50 — lOOK depending on the carrier density. 



We note that Pph{T) due to phonon scattering is always 
monotonic in T, and therefore the 'insulating' temper- 
ature dependence of p{T) around the low-density Dirac 
point cannot arise from phonon scattering which would 
always produce a metallic pph (T) increasing with increas- 
ing T (linearly at higher temperatures). It is, therefore, 
reasonable to conclude, as we do in this work, that the 
low-T insulating behavior of graphene p{T) around the 
low-density Dirac point, as observed experimentally in 
ref. [l3], is a result of the high-temperature and low- 
density (i.e. T/Tp < 0.5) screened impurity scattering 
phenomenon discussed in this work. More work will, 
however, be needed to understand this Dirac point in- 
sulating behavior quantitatively since the electron-hole 
puddle induced density inhomogeneity becomes impor- 
tant around the Dirac point. 
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APPENDIX A 

We provide, for the sake of convenience, in this ap- 
pendix (1) the various formula comparing the basic elec- 
tronic properties [e.g. Fermi wave vector, Thomas- 
Fermi (TF) wave vector, r^-parameter, density of states 
(DOS), Fermi energy, cyclotron frequency] in graphene 
and parabolic 2D systems; and (2) the analytic equa- 
tions for the electronic polarizability and (3) the resistiv- 
ity (arising from screened charged impurity scattering) in 
the low (T < Tp) and high (T > Tf) temperature hmits 
for graphene and 2D systems. The results are given in 
terms of carrier density (n), Fermi velocity vp (graphene) 
or effective mass m (2D systems), ground state degener- 
acy g (= 5s5u for spin and valley), background dielectric 
constant (k), electron charge (e), and Planck constant 
(h) in tables I, II, and III. 
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TABLE I: Electronic quantities. Note: The carrier effective mass (m) in 2D systems and the graphene Fermi velocity (vf) 
are assumed constant independent of carrier density (n) and defining the basic single-particle energy dispersion at wave vector 
q: e(q) = /ri'F|q| (graplicnc) or fi'^q^/2rn (2D systems). The degeneracy factor g = gsgv carriers the usual spin degeneracy 
(g^ = 2) and a valley (Ic^ciit^racy (f/,. = 2 fcjr graphene). 



(Quantity 



Parabolic 2D system 



Graphene 



Feriru wave vector (kp) 
TF screening wave vector (qxp) 
Wigner-Seitz radius (rs) 
DOS {D{E)) 
DOS at Ef (Do = D{Ef)) 

Fermi energy (Ef) 
Cyclotron frequency (wc) 
Landau level energy (Ei) 

Plasma frequency {LOp{q)) 



^/4Trn/g 
gme^ / kK^ 

gm/2%T{^ 
g'm/2'KfP 
2nnh^ /gm 
eB/mc 
{l+^)hL u, 1= 0,1,2, .. 



ge^ kp / kHvf 

e^/KhvF 
gE/{2'Kh^v%) 
gkp /{2'KhvF) 
hvF y/^im/g 
VF\/2ehBn 
sgn{l)vF^/2ehB\l\ 1 = 0, ±1, ±2, . 



TABLE II: Temperature dependent polarizability Yl{q,T). Note Do is DOS at Fermi energy. Here C^{x) is the Riemann's zeta 
function. 



Temperature (T) 



wave vector (q) 



2D n/Do 



graphene H/Do 

6 I Tp 



Low T ^Tf 



High T > Tf 



q = Q 
q = 2kF 

All q 



1 ^ (1 - v^)C(i)y^ 



i + /i(i-^)C(|)(^ 



t.l,--;- T 



T 



ln4 - 



, 2 



3/2 



TABLE III: Temperature dependent conductivity <j(T). Note ao = cr(T = 0). 



Temperature (T) 


2D system 


Graphene 


Low T <^Tf 
High T > Tf 


cT{T)/ao = 1 - 0{T/Tf) 
cr ~ 0{T/Tf) 


a{T)/ao = 1 - 0[{T/TFr] 
a - 0[(T/TFf] 



